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Abstract 



Qh! We study further the role of the boundary operator Ob for macroscopic loop 

r-| ! length in the stable definition of 2D quantum gravity provided by the [P, Q] = Q 

formulation. The KdV flows are supplemented by an additional flow with respect 
to the boundary cosmological constant a. We numerically study these flows for 
the m = 1, 2 and 3 models, solving for the string susceptibility in the presence 
5^ I of Ob for arbitrary coupling a. The spectrum of the Hamiltonian of the loop 

quantum mechanics is continuous and bounded from below by a. For large positive 
cr, the theory is dominated by the 'universal' m = topological phase present only 
in the [P, Q] = Q formulation. For large negative a, the non-perturbative physics 
approaches that of the [P, Q] = 1 definition, although there is no path to the unstable 
solutions of the [P, Q] = I m-even models. 
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In certain toy models of closed string theory, the [P, Q] = Q formulation[l] 
supplies the most general string equation compatible with the KdV flows: 

pu^^hjzn +\{nf = Q. (1) 

Here IZ = Z]fcLo(^ + \)tkTlk[p\: where the IZk are the Gel'fand-Dikii differential 
polynomials in p. A prime denotes vdjdz where v is the string coupling. The 
quantity z = —to and the normalisation chosen is TZq = 2, which fixes the other TZk 
by virtue of the recursion relation TZk_^.i = jT^k ~ \p ^fc ^ P^k ^^^ the requirement 
that they vanish when p = 0. 

Setting tk ~ 5fcm selects the string equation of the (2m — 1, 2) conformal minimal 
model coupled to 2D quantum gravity. In particular m = 2 is pure 2D quantum 
gravity. In this case the parameter z is the bulk 'cosnio logical constant', the chemical 
potential for the area of the random surfaces. In operator language, z = —to couples 
to the 'puncture operator', Oq in the action. Oq is the gravitationally dressed 
identity operator of the (3, 2) model and may be thought of as marking out a point 
on the surface. The other tk couple to other operators Ok in the theory. The 
function p{z) is the 'string susceptibility', the two-point function of Oq: p —< 
OqOq >= — r . r is the free energy of the theory. The presence of the KdV flows in 
the {2m — 1,2) models is most naturally phrased in terms of p: dt^p=< OkOoOo >— 
TZk+i- This is a statement about the organisation of the operator structure of the 
theory and flow between all the minimal models. Using this relation, the string 
equation and the recursion relation, all correlators of the local operators Ok may 
be determined. 

The [P, Q] = 1 definition of these models selects the solutions of (1) for the 
string susceptibility which satisfy 7^ = everywhere. Asymptotically, for the mth 
critical point, we have p ~ z^'"^ for z -^ ±oo. This leads to the familiar cases where 
for TO-even, real solutions are problematic [2] and for m odd, we have the BMP-type 
solution[3]. The [P, Q] = Q definition of these models selects the unique[4] solution 
of (1) with the asymptotics p ~ ^i/"^ for z -^ +oo and p ~ for z -^ — oo. These 
solutions provide a definition for p which is free from non-perturbative instabilities. 
Note that for m odd that the [P, Q]=Q solution for p differs from the corresponding 
[P, Q] = I solution, although they are both unique and well-defined. 
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Another distinguishing feature of the [P, Q] = Q definition is the existence 
of an TO = critical point. The [P, Q] = 1 definition does not possess such a 
point[5]. Setting tk = 0, {k ^ 0) in (1) yields the simple solution p = — z/^/(4z^). 
This solution for p does not support finite area surfaces and leads to recurrence 
relations which determine all operator products. Thus this m = critical point is 
a purely topological theory, in addition to the usual m = 1 topological point. This 
distinguishing feature feeds into all of the other rn-critical models in the [P, Q]=Q 
definition by supplying them with an interpretation in the 'weak coupling[5]' phase, 
the z —^ —oo limit. In this limit, the leading correction to the p = behaviour for 
any to is p = — z/^/(42^), and so the universal to = topological phase dominates 
the physics in this regime. 

The models may be studied further by the inclusion of macroscopic loops in 
the theory [6]. The features of the solution for p receive further attention here, 
for the discussion of macroscopic loop dynamics in these models reduces to ID 
quantum mechanics with Hamiltonian Ti = Q = u^ d^ — p. Here, the usual hermitian 
matrix model inspired prescription for the expectation value of the macroscopic loop 
wavef unction, 

< w{i) >= r < zle^^""^^^- - ^) |2 > dz (2) 



is adopted. For the [P, Q] = 1 definition, H has a discrete spectrum for the m~ 
even models. A simple exponential behaviour for < w{i) >, e~^o^ is obtained only 
in the large i limit, where eo is the lowest eigenvalue of 7i. Meanwhile for the 
TO-odd solutions the loop expectation diverges as £ ^ oo due to the presence of 
the exponential tail of the scaled charge density in the associated Dyson-gas-on-H 
problem[3]. Here again, the [P, Q] = Q definition distinguishes itself. The TO-critical 
models all have the same generic behaviour, in contrast to the [P, Q] = l definition. 
In considering the loop structure of the [P, Q]=Q theory, very simple arguments [7] 
yield the exponentially decreasing behaviour for loops. This arises naturally from 
the fact that the [P, Q] = Q definition has a realisation in terms of a Dyson gas on 
IEl_i_. This will be reviewed below. 

The length of loops is discussed most naturally by the introduction of the 
boundary length operator Ob to measure it. The associated coupling will be 
denoted by a, the boundary cosmological constant. The operator Ob in the 
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[P, Q] = l definition of these models was studied in ref. [8] . Its signature is identified 
by noting tfiat a may be set to zero by a redefinition of tfie couplings tk of the bulk 
operators Ok {k < m) in the mth model[8] and a shift in p. Perturbatively in the 
tfc, a is identified with the combination tm-i/tm for the mth model. Meanwhile 
Ob is a mixture of Om-i with all the lower operators 

m -1 

OB=Y.ik+0^trnOm-l+--- (3) 

k=l ^ 

where the dots represent higher order terms in tm-i- In the [P, Q] = Q definition 
the boundary operator and its coupling a arise in a more natural fashion than in 
the previous study. In ref. [7], the most general string equation in the presence of 
the coupling a and compatible with the KdV flows was derived: 

(p-a)n'-^nn" + ^{nf = o (4) 

which, by virtue of the fact that the first derivative of (4) is a scaling equation, 
yields 

The equations (4) and (5) arise naturally in the double scaled limit of a Dyson gas 
on 1R_|_ [7] [9] . The parameter a represents the scaled position of the infinite potential 
'wall' defining the 1R+ topology. 

The identification of Ob is most easily made[7][9] by observing that equations 
(4) and (5) are invariant under the infinitesimal transformation: 



p- 


^ p = p-e 


a - 


-^ a = a — e 


z - 


^ z = z - e-ti 


tk- 


3 

^tk=tk + e{k + -)tk+i 



(6) 



which is a Galilean transformation, supplemented by a transformation for a. The 
redundancy of a is manifest here, as it may be set to zero by performing a Gnite 
Galilean transformation, which redefines all the couplings tk and shifts p by —a. At 
the same time, using equation (2) for the expectation value of the macroscopic loop 
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wavefunction, we see that the resulting shift in p displays the e"*^ loop behaviour. 
This completes the identification of a as the boundary length operator's coupling. 
Meanwhile equation (5) yields the precise combination of the Ok which make up 
the boundary operator: 



oa Ja da 



dz 
n ' da 

< w{e) > + E (^ + 2^^^drT < "'^^^ > 

k=l "- ^ 



from which we deduce: 



— <w{i)>=<OBw{i)>=i<w{i)> (8) 

where 

is the insertion of the boundary operator. We thus have the Ward identity 

kObU w{e,) >= E ^j < n ^i^^) > (10) 

i j i 

which is usually referred to as the L_i Ward identity. Indeed in the language of 
Virasoro constraints, equations (5) and (4) may be written as the familiar L_i 
and Lq constraints, with a boundary term due to the presence of a: L_it = d^^r; 
Lqt = adcfT. Similarly, the rest of Virasoro constraints appear as the familiar ones 
modified by boundary terms: LnT = a^^^d^T. So in the [P, Q] = Q definition, the 
redundancy of the coupling a turns into a redundancy of the L_i constraint, as the 
finite Galilean transformation above amounts to the following automorphism of the 
constraint algebra[7][10]: L„ =e~'^ ^L^e^ i, under which the L_i constraint 
disappears, and we are left with L^t = (n > 0). 

In ref.[9], the loop equations for the [P, Q] = Q definition of the (2m — 1, 2) 
models were derived and studied. Here, we turn our attention to the study of 
equations (4) and (5). For a = we have the previously studied cases: there is a 
unique pole-free solution^ for the mth. model with the asymptotics p ~ ^i/"^ for 



^ This has been demonstrated for m=l,2 and 3, and conjectured to be true for 
TO > 3. 



z -^ +00 and p = for z -^ —00. These solutions have been studied analyticaUy 
and numericaUy elsewhere[l][4][5]. For non-zero a, an asymptotic analysis of the 
string equation is very similar to the case for cr = 0. It will not be repeated here. 
The leading exponential corrections to the asymptotics p = z^'"^ for z -^ +00 and 
p = (J for z -^ —00 can be calculated using a WKB ansatz. The coefficients of these 
exponentials are all determined by the boundary conditions. Therefore all of the 
integration constants have been determined locally, and there is at most a discrete 
number of solutions with the above asymptotics. A consequence of this is that the 
only infinitesimal change that we can make to the solution is via KdV fiows in the 
tk or the cr-fiow with equation (5). So for any rn-critical model, local solution space 
is spanned by one-parameter KdV fiows to other models, and the fiow in the single 
parameter a given by equation (5). Therefore, beginning at any of these unique 
solutions, any solution obtained using these fiows is also unique. 

The KdV fiows between models were studied in ref.[4], where numerical meth- 
ods were used to display intermediate flows, the solutions to the interpolating string 
equation (1). Here, the accompanying a-fiows in the string equation (4) are studied. 
The solutions were found primarily using the FORTRAN NAG finite element library 
routine D02RAF. This routine allows the solution of several simultaneous first 
order ordinary differential equations with given boundary conditions by use of 
Newton interpolation from an initial approximate solution supplied by the user. The 
equation ^TZ'" — \plZ — {p — a)lZ = was studied, which is the first differential 
of (4) with a factor TZ removed. This enabled better calculation of the numerical 
derivatives: solving (4) directly would involve the implicit calculation of 1/7^, which 
diverges for large positive z, as 7^ = in that asymptotic limit. The above equation, 
a differential equation of order 2m + 1, was transformed to 2?ti + 1 first order 
equations. Boundary conditions were taken from the above asymptotics. 

The first case studied in this way was m = 1. The equation was solved with 
TZ = ^p — z on a mesh of 2000 points. The routine converged rapidly to the solution 
for each value of a with an error of better than 10~^. The results are shown in 
figure 1. The solution for m = 2 with a positive was carried out in the same way. 
Here TZ = f ( — |p + P^) — -z, and a mesh of 2000 points was used. The results are 
displayed in figure 2. For a negative, however, the solution becomes progressively 
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more oscillatory, and the derivatives take on increasingly large values. Eventually 
the required accuracy could no longer be obtained using D02RAF. Here, the a~ 
flow equation was employed to make further progress. In order to solve for more 
negative a, the solution for some previous value of a is used as a starting solution. A 
numerical integration of equation(5) is then carried out using Euler's method. This 
gives the next initial approximation for D02RAF, and the process can be iterated to 
give solutions for more negative a. The results are shown in figure 3. Unfortunately, 
at a = —1.25, both rounding errors and systematic errors become so large (of order 1, 
and much larger for higher derivatives) as to prevent convergence. Rounding errors, 
and the reliance on high derivatives, prevent the use of more accurate integration 
methods such as Runge-Kutta. However, the pattern of behaviour of the solution 
for increasingly negative a has been established. The m = 3 solution was found in 
a straightforward way, using the method described for the m = 2, a negative case 
above. Here however, there were no problems with convergence, and the program 
required relatively few mesh points to give small errors (typically 2000 for error of 
< 10~^, against 15,000 for the m = 2 case for a < 0). Results are shown in figure 4. 

When studying these curves, the associated double scaled Dyson gas problem 
should be recalled[ll][12][l][7][9]. In scaled coordinates A^, in the spherical approxi- 
mation, the charges are concentrated on a single cut extending along the semi-open 
interval [p, oo). There is an infinite potential wall at position As = a, restricting the 
problem to 1R+. 

Without the wall, the topology is IR which leads to the critical physics of the 
[P, Q] = l definition. The critical physics of the model is entirely characterised by 
the local behaviour of the charge density in the neighbourhood of its endpoint[13], 
at position As = p. The density of charges vanishes at least as a square root at 
its endpoint, and for the mih. critical model it has to — 1 extra zeros there. The 
resulting string equation derived for p is 7^ = 0. In examining the stability of 
the critical models a study at the spherical level of the effective potential for one 
eigenvalue reveals that for an odd number of zeros it is energetically favourable for 
eigenvalues to leave the single cut and move to a new configuration[9][l]. This is 
the origin of the m-even models' instability. 

In the [P, Q] = Q formulation the effective potential described above is sup- 
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plemented by the infinite potential wall, which has the effect of stabilising the 
probleni[l]. In the asymptotic analysis at large 2, there is the following behaviour: 
For z -^ +00 the density pulls away from the wall, and the neighbourhood of the 
endpoint is the same as that of the [P, Q] = 1 definition. Therefore the genus by 
genus physics for large positive z is identical to that of the [P, Q\ = \ definition. At 
the string equation level, this is realised by the fact that the solution of equation 
(4) satisfies 7?. = in this limit. For ^ ^ — 00 the density pushes up against the wall, 
and the square root zero is replaced by an integrable square root divergence. This 
leads to p = cr for the leading behaviour of the endpoint. 

The non-perturbative positions of the charges in the Dyson gas are directly 
related to the spectrum of the hamiltonian Ti = Q = f^ d^ — p. For the zero a case, 
the single cut nature of the charge density is consistent with a continuous spectrum 
for Ti with no discrete states below the continuum. This may be confirmed by 
examining, for example, the form of the solution for p for m = 1. The small well is 
too shallow to support such discrete states. This is immediately true for all of the 
TO-critical solutions for p: they are connected by global KdV fiows[4] and the KdV 
fiows are the isospectral deformations of the hamiltonian 7i. 

Moving on to discuss the non-zero a case, we first consider the nature of 
the cr-fiows. They represent the first in a series of important non-isospectral 
deformations of 7i, the Galilean transformations. The next in this series are the 
scalings. Demanding invariance under these transformations leads to the L_i and 
Lq Virasoro constraints. The rest of the Virasoro constraints follow from higher, 
non-local transformations which may be derived using the recursion operator for 
the KdV hierarchy. Using Galilean transformations, the continuous spectrum of 7i 
cannot be deformed to a spectrum with any discrete states. This is consistent with 
the Dyson gas picture where the Virasoro constraints arise from diffeomorphisms 
in IR_|_, with the terms due to a arising from varying the position of the wall. By 
definition, diffeomorphisms preserve the single cut nature of the charge density, and 
thus the bounded continuous spectrum of 7i. In particular, L_i which is the cr-fiow 
at the r-function level, is a translation of the Dyson gas and hence just a shift of 
the spectrum of H. 

The case of to = 1 is particularly simple. Starting with a = 0, the Galilean 
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transformation (6) may be used to move to non-zero a. The parameter ti remains 
invariant, and it may be set to 1 without loss of generahty. The result of this process 
at some finite a (positive or negative) is a shift in p and z by a. This corresponds 
to translating the whole cr = curve along the diagonal line p = z. This is the 
process illustrated in figure 1, which is the solution of (4) at m = 1 for arbitrary 
values of a. Away from 171 = 1, the processes involved are not as simple. This is 
because the Galilean transformation inevitably switches on couplings tk {k < m) to 
operators other than the puncture operator Oq, Om and Ob- For a going positive 
at fixed z the wall moves to the right, in scaled coordinates As, and pushes up 
towards the density. This is the reverse of the z ^ —oo limit for fixed a discussed 
before and asymptotically p = a — u'^/{4z'^) + . . .. In this large a regime, surfaces 
of finite area are increasingly suppressed and the model is dominated by the 171 = 
topological phase. See figures 1 2 and 4. For a increasingly negative, the wall 
moves to the left, and the solution for p gradually begins to resemble the [P, Q] = 1 
?7i-critical solution, although this process can never be completed for finite a. In 
particular, for the m = 2 (pure gravity) case, the solution develops progressively 
sharper oscillations for negative sigma. Naively, this suggests that in the limit 
a — > — oo, the oscillations evolve into the poles of the Painleve I equation. However, 
this is an ill-defined limit, as can be seen from many points of view. There is no 
way to continuously deform to such a solution of Painleve I in the same sense that it 
is not possible to connect the BMP solution for the [P, Q] = l definition of m = 3 to 
it via the KdV fiows[12]. Furthermore the spectrum of Ti for the singular solutions 
of Painleve I is discrete, and there is no way to deform the continuous spectrum 
of H for solutions of (4) into a discrete one via cr-fiow, as discussed above. This 
generalises to all of the m-even cases. For m-odd, the spectrum is continuous and 
not bounded from below for the BMP-type solutions of the [P, Q] = 1 definition, 
and it is conceivable that they are connected to the [P, Q] = Q m-odd solutions 
in the a -^ — oo limit. This picture is possibly consistent with one in which the 
wall has moved away sufficiently far to the left so as to have a vanishing effect 
on the non-perturbative physics. For the m-even models, the reduced presence of 
the wall in this limit allows the non-perturbative physics to more resemble that of 
the [P, Q] = 1 definition, but it can never be completely removed so as to recover 
instability. 
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Figure Captions 

Fig. 1: The cr-flow in the m = 1 model for positive and negative a. 

Fig. 2: The cr-fiow in the m = 2 model (pure gravity) for positive a. 

Fig. 3: The cr-fiow in the m = 2 model (pure gravity) for negative a. 

Fig. 4: The cr-fiow in the to = 3 model for positive and negative a. 
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